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The radiative decay is believed to be an ideal lab to study hadronic structure of newly observed 
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I. INTRODUCTION 



In the field of heavy hadrons, it is noticed that some experimental observations obviously 
deviate from our theoretical predictions, so definitely, such "anomalies" need to be clarified. 
Supposing the experimental measurements are right, there must be some loopholes in our 
present understanding of the nature, either there exists a contribution from new physics 
beyond the standard model (SM), or the concerned hadrons have exotic structure such as 
hybrid, multi-quark etc. It is believed that radiative decays on this aspect provide an ideal 
lab to testify the hadron structures and as well help to identify the quantum numbers for 
a newly observed resonances. Even though, the reaction mechanism for radiative decays is 
free of strong interaction (at tree level), the non-perturbative QCD effects are still involved 
in the hadronization of emerging hadrons, as well as the wavefunction of the parent hadron. 
To search for new physics or possible exotic components of the hadron, one must thoroughly 
study how the non-perturbative QCD affects the decay rates in a reasonable theoretical 
framework. In this work, we are going to derive the formulas for the radiative decays of ++ 
and l +_ heavy mesons which are supposed to have regular qq' structures, in the light front 
quark model (LFQM). Through a comparison of the theoretical predictions with data, the 
consistency degree would reveal if the concerned hadron possesses the regular qq' structure 
or has an exotic component. 

Duringrecent years several new particles such as D s (2317)[l ; ], D s (2460)0, 3, X(3872)0, 
X(3940)ja| and y(394O)0] have been observed in experiments and re-confirmed as new 
resonances. However it is not easy to identify their inner structures, i.e. if the constituents 
of the resonances are indeed just q and q' or with something else (gluon or qq pair etc.). 
Some of them are speculated as exotic states. The reasons may be twofold: one is their 
peculiar decay modes whereas the other is that the theoretical expectations on their excited 
states are not well consistent with data or even missing. For example when D s (2317) was 
announced it was considered as a molecular state or a tetraquark, however it may just be 
a regular p— wave cs meson with J p = + as it appears in particle data book0. Because 
radiative decay is fully governed by the electromagnetic interaction, the reaction mechanism 
is relatively simple and the emitted photon can be well measured in experiment, it may offer 
a good opportunity to justify the quantum numbers and constituent structure of the newly 
observed resonance. Namely, one may check whether with the simple qq' assignment the 
theoretical prediction agrees with data. 

For p— wave particles there are three degenerate states ++ , 1 ++ and 2 ++ with the total 
intrinsic spin S = 1 and one singlet l + ~ with S = 0. It is well known that Xco, Xci an d Xc2 
and h c are triplets and singlet p— wave charmonia respectively. 

In this work we will calculate the rates of ++ —> 1 7 and l + ~ — > CT + 7 in LFQM 0-0 



which has been successfully applied to evaluate rates of semileptonic and non-leptonic decays 
of s— wave heavy mesons. For the p— wave mesons, the wave functions of ++ , 1 , 1 H and 
CT + have been constructed and the leading Feynman diagrams are simple, with them we are 
able to calculate the corresponding transition amplitudes. 

Since the branching ratio of Xco J 7^7 is wen measured, we first calculate B(x c o ~^ 
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J/ip^i) in LFQM to fix the model parameters and check the validity degree of this approach, 
then with the formulas for the transition ++ — > 1 7 we estimate the rates of xw T(1S')7 
and .D s (2317) — > D* + 7. These formulas can also be applied to study decays of other ++ 
states. Recently the spin singlets h c and hb attract intensive interests of both experimentalists 



and theorists [16Ml8j . By the formulas for l + ~ — > 0~ + 7 we calculate the widths of h c — > 7/ c 7 
and hb —> rjbj- The results can be compared with the data which will be available soon 
at the BES II and future B-factory. This comparison definitely assures us if they are pure 
p— wave heavy quarkonia. 

After the introduction we derive the formulas for the transition ++ — > 1 7 and l +_ — > 
CT + 7 in section II. Then in Sec. Ill, we numerically evaluate the decay widths of Xco ~ ► J/^li 
Xto Y(1S)7, D s (2317) — > D* + 7, h c — > r]^ and h b —> rj b ^ and make some discussions. In 
the last section we give a brief summary and discussion. Some notations and definitions of 
relevant quantities are collected in the attached appendix. 



II. THE FORMULA FOR THE DECAYS 0++ -> 1 7 AND 1 



(r+7 



In Ref.pjj] Chung studied various transitions and derived the corresponding amplitudes, 
for ex amp le, he determined the amplitude structure of bi (1235) — > coir. For transition l +_ — > 
_+ 7 19j a photon exists in the final states, thus the gauge invariance demands the transition 
amplitude to be in the form 



A = F * (g 



w - k . q ) £ > 



(1) 



where q and k represent the momenta of the photon and daughter meson, e' and e are 
the polarizations of the vector-meson and photon respectively, and F is the form factor 
which is what we are going to derive and numerically compute in this work. Though the 
++ — > 1 1 was not discussed in Ref. 19], its amplitude structure is the same as Eq.([T]) 
and it can be seen from analyzing the J p characters of the involved mesons. 

Below, following the schemes given in literature, we will calculate the transition ampli- 
tudes of interest in LFQM. 



A. the decay of ++ — > 1 7 



The vertex functions of ++ and 1 are respectively 12 



-Ms 



(2) 
(3) 



where Hs, Hy and Wy are defined in Ref. 12| and p a and pb are the momenta of constituents 
of the corresponding meson. 
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FIG. 1: Feynman diagrams for radiative decay. 



The transition amplitude corresponding to the left diagram of Fig.l is written as 



A a = ieei 



(2tt)< 



4„ _HsHv_ a lu a* 

NiNiNz ^ 



(4) 



where 



= Tr{(-^ 2 + m 2 )[ju - Pl Tzr P2 " KA + m i)lti{^i + m i)}, 



in/ 



P 



i2 



Wy 

m\ + ie, N 2 = p 2 — m\ + ie, eei( 2 ) is the electric charge 



JVi = - m 2 + ^ jy/ 
of the quark of u- or d-types. In the light front frame, Pi is decomposed as (p^,pt,Pi 
Integrating out Pi with the methods given in Ref. U] one has 



with 



,4 H S H V IU r 2 h s h v Aa - 1/ A „ 
a Pi — r— — s" e -> -ct / dxid p± ; — —s,„,e' 



/i s = (M 2 -M 2 )< 



1 M 2 



/i v = (M' 2 - Af'o 



iV c V2M 2V3M (3 
1 



I x'ix'2 



N- 



(') 



J' 



(j>'(nS), 



(5) 



where M and M' represent the masses of decaying and produced mesons and the relation 



between e(e') and can be found in the appendix of Ref. 12 



To include the contributions from the zero mode, p ifl , Pu L ,Pi fJi Pi u and Wy in s a should 
be replaced by appropriate expressions as discussed in Ref.|l2|], for example 



Wy — > Wy = M + Uli + m 2 

_^ xi .x p±q± 

Pin -> + ^2~ ~tf~> q ^ 



(6) 
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with V = p + k and p is the momentum of decaying meson. 

The definitions of Mo and Mq are presented in the appendix, as more details about the 
derivations and relevant notations can be found in Ref. 12j. In this framework, s" replaces 



s" and is written into the form 



kpSlv 
k ■ q 



(7) 



with 



-2 



mi 



(M 2 - M' 2 - #1 + iV() + m 2 (Nx +N[-q : 



4 A12 (2 mi 2 + Ami m 2 + 2 m 2 2 — M 2 — M' 2 + g 2 — 2 m-i — 2 m 2 iiy + 2Z2) 



4 fc • q A { 2 ] (-Ni - N[ + q 2 ) 4 k ■ q (Ni + N{ - q 2 + m x w v 
F 2 = 1 h 

w v w v 
Ak-qA { 2 } {-2 [2 (mi + m 2 f - M 2 - M' 2 + g 2 ] +4mi% + 4m 2 %} 



4/c-g^ 2) {-2 [2 (mi + m 2 f - M 2 - M' 2 + g 2 ] +4m 1 tOi/ + 4m 2 %} 



+ 



+ 



4 jfe • g aS 1} 4 (mi + m 2 ) 2 - 2 M 2 - 2 M /: " - N x - N[ + 3 g 2 - 2 (3 m l + m 2 ) w v 



/2 



and A^\i = 1 ~ 4, j = 1 ~ 4) are presented in the attached appendix. 
We define the form factors as following 



/i(mi,m 2 ,ei) 



/ 2 (mi,m 2 ,ei' 



eei /• , 2 Fi00'M o 
ax 2 ct 



32tt 



3^ 1 dX2£p± 



F 2 (jxj)'M 
V6(30M o M' o x 1 ' 



(9) 



where and <p' are respectively the wavefunctions of the initial and final mesons. These form 
factors will be numerically evaluated in next section. 

With these form factors the amplitude corresponding to the left diagram of Fig.l is 
obtained as 



A a = fi(m 1 ,m 2 ,e 1 )e- e' + f 2 {m 1 ,m 2 ,ei 



k ■ e'q ■ e 
k ■ q 



(10) 



The right diagram is just the charge conjugation of the left one of Fig.l, so that one can 
immediately write it down 



A b = f 1 (m 2 ,m 1 ,e 2 )e ■ e' + f 2 (m 2 ,m 1 ,e 2 



k ■ e'q ■ e 

t 

k • q 
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The total amplitude is simply a sum of A a and A b : 

A =A a + A b . (12) 

Comparing Eq. (1) with Eq. (12) we determine the full form factor as 

F = /i(m 1 ,m 2 ,ei) + fi(m 2 ,m 1 ,e 2 ) = -[/ 2 (mi, m 2 , e x ) + / 2 (m 2 , mi, e 2 )]. (13) 

In principle one can calculate fi = /i(mi, m 2 , ei) + /i(m 2 , m 1 ,e 2 ) and / 2 = / 2 (mi, m 2 , ei) + 
f 2 (m 2 , mi, e 2 ) separately and fix the form factor, then go on obtaining the decay rate. 



B. the decay rate of 1" 



(r+7 



The vertex functions of + and 1 + are presented in Ref. 



12 



they are respectively 



-iHi 



1 A 



(14) 
(15) 



where Hp, Hi a and Wi A are some relevant functions which are slightly lengthy and can be 
found in Ref.[12j, for saving space, we do not repeat them here. 

The transition amplitude corresponding to the left diagram of Fig.l is written as 



A a = eei 



Nr. 



where 



[IV 



Tr{(-^ 2 + m 2 )[ 



(P'l ~V2)v 



„ Hp Hi ^ , 



(16) 



W: 



1 A 



7bWi + mihufa + m i)7s}- 



Integrating out p 1 we have 



4 H P Hi A IV . f 2 /ip/li A 



where 



NiN[N 2 



hi. 



hi 



(M 2 - m 2 ; 



(M /2 -M' 2 )< 



'^iX 2 



x 2 N x N{ 
y/2 



s a e' 



(17) 



1 



0(nS), 



I X\X 2 



^</>'(nS). 



°1 N c ^2M' 

With a replacement similar to Eq.(j6]), s* is re- written into the form of Eq. ([7]) with 

4g 2 8Z 2 . 



Ft 
F 2 



(2) 



16 mi m 2 8 m 2 2 4 M 2 4 M' 2 



2£;-g[4 : 



,27V, 27V{ 



2g\ 



U>i A Hiiyi U>l A 101^ 

2 Nt 2N[ 2q 2 



A 



Wi A Wi A Wia 

( 2 ) 8 mi 2 16mim 2 8m 2 2 



+ 



4 X) v 

~ Wiyl 101^ Wia' 

4 M 2 4 M' 2 " ~ 2 



WA 



Wi A 



A 



(2)/— 8 mi 2 lQrriim 2 8m 2 2 4M 



+ 



4M' 



4g 2 

Wia' 



/2 



4^ 

Wi A ' 



)]■ 
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TABLE I: the form factors f\ and /2 



decay mode 


h 


h 


width (keV) 


Bft(the) 


B^(exp) 


Xco -> ^7 


-2.83 ±1.37 


0.89±0.09 


152±31 


(1.46±0.31)% 


(1.17 ±0.08)% [7] 


Xw -> T(15)7 


-0.50±3.01 


0.82±0.08 


21.3±4.7 




(1.76±0.35)%[7] 


D s (2317) -> D s * 7 


-1.11±0.22 


0.25±0.03 


17.1±3.9 


> (0.45 ±0.11)% 




/i c ->■ ??c7 


-2.36±0.34 


2.47±0.18 


685±89 


> (69 ± 9)% 


(51±6)%[7] 


^6 ->■ %7 


-2.40±0.31 


1.78±0.20 


36.9±8.7 




(49.2±5.7i|J)%[20] 



TABLE II: Predictions made in terms of various approaches and experimental data (if available) 
(in units of keV) 



decay mode 


[24] 


[25] 


our results 


exp 


XcO -> J 1^1 


121 


input 


152±31 


122±12 [7] 


Xbo -> T(15)7 


29.9 


85±4 


21.3±4.7 




-> »7c7 


560 


634±32 


685±89 




/ib -> r?67 


52.6 


271±14 


36.9±8.7 





III. APPLYING THE APPROACH TO ANALYZE RADIATIVE DECAYS 



In this section we first test the formula by comparing our theoretical evaluation on the rate 
of a well measured decay mode with the data and confirm its validity, then apply it to predict 
the rates of radiative decays of ++ and 1 H which will be experimentally measured soon. 
We select transition Xco — > J / ^7 as a probe to check the approach since its branching ratio is 



well measured. Setting m c = 1.4GeV and the model parameters /3^( Xc ) = 0.631GeV|21 



we 



get the form factors j\ and ji for the transition Xco ~^ J/^l which are presented in Tab. I. 
One can immediately notice that fi is sensitive to the variation of the parameters, especially 

m c , but /2 is insensitive. The problem originates from the fact that iVf is proportional 

to a cancelation of two large numbers. Even though Nj resides in both F\ and F 2 , its 
contribution is suppressed by wy or wi^ in F 2 , so that F 2 is insensitive to the parameters. 
For reducing uncertainties of our theoretical computations, we will calculate F in terms of the 
relation between f 2 and F in Eq. ffT3l) . Our theoretical estimate of the width of Xco J /''Pi 
is (152±31)keV and the corresponding branching ratio is (1.46±0.31)% which is consistent 
with the data (1.17 ± 0.08)%|7j within a tolerable error range. 

Then we study the transition Xbo — > T(lS)'~y using the parameter m& = 4.8GeV and 
(3y( Xb ) = 1.288GeV which were fixed in Ref. [221] by fitting other well measured channels. Our 
estimate of T{ X bo -> T(1S)7) is (21.3±4.7)keV. With the branching ratio (1.76±0.35)%f" 
the total width of Xbo is estimated to be (1.21±0.36) MeV. 

Supposing _D S (2317) to be a regular bound state of cs and setting m s = 0.37GeV, (3d s 



We vary the parameters within a ±10% range to estimate the corresponding errors. 
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0.592GeV[23|, we can calculate the rate of -0^(2317) — > D* + 7 using the above formulas. 
We obtain its partial width as (17.1±3.9)keV. Namely, if -D s (2317) is a ++ regular meson, 
possessing the p— wave cs structure, our numerical prediction should be consistent with the 
data which will be available at BES II or B-factory soon, therefore the consistency degree 
with data can help to confirm or negate its assignment cs. 

In fact these formulas can be applied to study radiative decays of other heavy mesons 
such D , B , B s0 and B c0 . 

Once h c and hb were experimentally measured, their special behaviors draw intensive 
interests of theorists. The main focus is if they are regular quarkonia or have exotic compo- 
nents. Using the formulas for l +_ — > CT + 7, we calculate the branching ratios of h c — > rj^ 
and hb — > 77^,7 and the numerical result are presented in Tab. I, which can be used to an- 
alyze their characters. In terms of the measured upper limit of the total width of h c we 
estimate the branching ratio of B(h c — > 7^7) > (69 ±9)% which is a little larger than the 
data (51 ± 6)%. If h c (hb) is a pure charmonium (bottomonium) our estimate indicates that 
its total width should be around 1 MeV (80keV). To further identify the inner structures of 
h c and hb, more precise experiments are needed. 

It is noted that the masses of initial and finial mesons we used in our numerical compu- 
tations are taken from Ref. j3] except the mass of rib which is 9402. 4MeV 20]. Some parallel 
researches on these decays can be found in Ref . 24|, |25J and a comparison of their results with 
ours is made and listed in Tab. [TTJ Because of large errors in the inputs for our numerical 
computations, our predictions are of relatively large uncertainties. To further testify the 
model or constrain the model parameters, much more precise experiments are needed. 



IV. SUMMARY 

In this work we derive the formulas for the radiative decay of ++ and l +_ heavy mesons 
and numerically compute the rates. We formulate the transition matrix elements and ex- 
tract the form factors in the LFQM. To check the validity degree of this approach where 
the model parameters were fixed by fitting other physical processes in previous works, we 
calculate the decay widths of Xco ~^ J/i>l an d compare it with data. Considering both 
theoretical and measurement uncertainties, the theoretical result on the branching ratio of 
Xco -> J/ipl (1.46±0.31)% is satisfactorily consistent with data (1.17 ± 0.08)%. With the 
same scenario, we predict the width of xw — > Y(lS)j to be about (21.3±4.7)keV and the 
width of _D S (2317) — > D* + 7 to be (17.1±3.9)keV. Comparing those results with the data 
which will be available at BES and B-factory, one can further confirm the inner structure of 
Xbo and £> s (23 17). 

In terms of the data B(xbo — > T(lS)"y) we estimate the total width of xto a s 
(1.21±0.36)MeV which would be easier to be experimentally checked. As for the l +_ de- 
cay, we evaluate the decay rates of h c — > r/ c 7 and h b — > 77^7 which are (685 ± 89)keV and 
(36.9 ± 8.7) keV respectively. These formulas deduced in this work can also be applied to 
study radiative decays of other ++ and l +_ particles. 
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Appendix A: Model description 



In the conventional light-front model, a meson containing a quark qi and an antiquark q 2 
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with its total momentum P and spin J can be expressed 

\M(P 2S+ \Lj,J z )) = f{d 3 p l }{d 3 p 2 }2(2nf5 3 {P-p 1 -p 2 ) 



x ^ls (Pi,P2, Ai, A 2 ) | qi(pi, Ai)g 2 (>2, A 2 )>, 

Ai,A2 



(Al) 



where ^/'[ J s z (pi,p 2 , Ax, A 2 ) is wave function in momentum-space, Ai and A 2 denote helicities, 
Pi, p 2 are the on- mass-shell light-front momenta defined by 



Pi = (ptiPi±)i Pi± = (PnPi)i Pi 



m 2 + p 2 ± 



P 



+ ' 



and 



{d 3 Pi } 



dpfd 2 pij 
2(2vr) 3 



S 3 (Pi) = S(pt)S 2 ( PlL ). 



(A2) 



(A3) 



From the eigen-equation 



26 



H LF \M(P 2S+1 ,Lj, J z )) = ^r^\M(P 2S+ \Lj,J z )), 



one can deduce a light-front Bethe-Salpeter equation 

dx'd 2 p'j_ 



(M 2 -M o 2 )0(x,p ± ) 



(x',p> ± ). 



(A4) 



(A5) 



In principle by solving Eq. (lA5j) one can obtain the momentum distribution amplitude 
<f)(x,p±). However the two-body interaction kernel V e // is complicate so some phenomeno- 
logical amplitudes are chosen in practicable theoretical calculations. In this work we use the 
Gaussian-type amplitudes, for example 



7T 



^) = 4(^ V(fa _ 



dpz ( P 2 z +p\- 
exp( 



(A6) 



where f3 is the model parameter to be fixed by calculating the decay constant. A parallel 
approach about LFQM was done in Refs. 27J, |28( where the Gaussian-type amplitude is 
adopted as a trial function for the variational computations of the QCD-motivated effective 
Hamiltonian. 
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Appendix B: Notations 



Here we list some variables and notations appearing in the context. The incoming meson 
in Fig. [TJhas momentum P = P1+P2 where p\ and P2 are the momenta of the off-shell quark 
and antiquark and 

p+ = Xl P + , p\ = x 2 P + , 

PlX = XiPj_ + p±, P2J_ = x 2 Pj_ - p±, (Bl) 

with Xi and pj_ are internal variables and x\ + x 2 — 1. 
The variables Mo, Mq and N± are defined as 



M 2 _p\ + m\ t p\ + ml 



Xi x 2 



M = a/Mo 2 - (mi - m 2 y. (B2) 



witn p z - 2 2x2Mo . 

Aij(i = 1 ~ 4, j = 1 ~ 4) are 



/.(!)_/. (1) P-L • /. (2) _ „2 (P± ■ <7±) 



A x ® = ^, A 2 ^ = A^-^^, 



2 



A 2 ^ = {A^)\ A^=A^A 2 ^\ A 4 W = (A 2 «) 2 - 
Ai (3) = Ai^, A 2 < 8 > = A 2 (1) Ai (2) , A 3 ® = Ai (1) A 2 (2) , 
A 4 ® = A 2 ^A 2 ^\A^ = (Alpl, a 2 ^ = A^A^\ 

A^ = A^A^\ A^ = A^A^ - (B3) 
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